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In  this  work,  density  functional  theory  (DFT)  calculations  have  been  used  to  investigate  chemical 
sensing  on  surfaces  of  single-layer  M0S2  and  graphene,  considering  the  adsorption  of  the  chemical 
compounds  triethylamine,  acetone,  tetrahydrofuran,  methanol,  2,4,6-trinitrotoluene,  o-nitrotoluene, 
o-dichlorobenzene,  and  1,5-dicholoropentane.  Physisorption  of  the  adsorbates  on  free-standing 
surfaces  was  analyzed  in  detail  for  optimized  material  structures,  considering  various  possible 
adsorption  sites.  Similar  adsorption  characteristics  for  the  two  surface  types  were  demonstrated,  where 
inclusion  of  a  correction  to  the  DFT  functional  for  London  dispersion  was  shown  to  be  important  to 
capture  interactions  at  the  interface  of  molecular  adsorbate  and  surface.  Charge  transfer  analyses  for 
adsorbed  free-standing  surfaces  generally  demonstrated  very  small  effects.  However,  charge  transfer 
upon  inclusion  of  the  underlying  Si02  substrate  rationalized  experimental  observations  for  some  of  the 
adsorbates  considered.  A  larger  intrinsic  response  for  the  electron-donor  triethylamine  adsorbed  on 
M0S2  as  compared  to  graphene  was  demonstrated,  which  may  assist  in  devising  chemical  sensors  for 
improved  sensitivity.  ©  2074  AIP  Publishing  LLC.  [http://dx.doi.org/10.1063/L4871687] 


I.  INTRODUCTION 

Two-dimensional  atomic  crystals  have  recently  emerged 
as  an  alternative  to  graphene,  showing  improved  attributes  in 
some  cases.*  Transition  metal  dichalcogenides  have,  in  partic¬ 
ular,  drawn  attention  for  application.^  Recently,  semiconduct¬ 
ing  single-layer  M0S2  (SL-M0S2)  field  effect  transistors 
(FETs)  demonstrated  some  improvement  in  ON/OFF  current 
ratio  and  relatively  high  mobilities,^”®  however,  transport  could 
be  complicated  by  localized  states^  or  atmospheric  adsorbates.* 
On  the  other  hand,  for  practical  applications,  M0S2  FETs 
emerged  as  promising  for  chemical  sensing,  e.g.,  of  NO,^’***  or 
NO2  and  NH3.**  Changes  in  conductance  upon  chemical 
adsorption  on  SL  M0S2  flakes  in  comparison  to  graphene  for 
various  analytes  demonstrated  a  robust  response.*^  A  relatively 
large  change  in  sensitivity  was  noted  for  the  electron-donor 
triethylamine  (TEA)  by  Perkins  et  al.}^  while  the  response  to 
electron-acceptors  was  minimal.  SL-M0S2  FETs  are  known  to 
exhibit  n-type  doping  although  the  origin  of  it  remains  mostly 
unclear,  presumably  due  to  impurities.**  Unlike  SL-M0S2, 
SL-graphene  is  known  to  exhibit  p-doping.  *"*’**  Indeed,  as  we 
previously  pointed  out,*®  understanding  the  behavior  of  a 
nano-FET  is  complicated  by  various  scatterers,  thus  necessitat¬ 
ing  understanding  of  the  intrinsic  mechanism  of  an  analyte’s 
adsorption,  which  can  achieved,  in  part,  by  computational  pre¬ 
diction.  Because  M0S2  could  become  a  two-dimensional  mate¬ 
rial  of  choice  for  chemical  sensing,  we  aim  in  this  work  to 
systematically  examine  parameters  that  may  affect  the  elec¬ 
tronic  structure  upon  molecular  adsorption,  both  on  M0S2  and 
for  comparison  on  graphene  surfaces.  Theoretical  work  has 
been  limited  so  far,  with  a  few  examples  for  either  adsor¬ 
bates*”*’*^  or  interaction  with  the  substrate*®’**”^**  for  one  type 
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of  surface,  while  M0S2  vs.  graphene  has  not  been  consistently 
assessed  for  chemical  sensing  application. 

In  this  work,  we  applied  density  functional  theory 
(DFT)  for  a  comprehensive  and  consistent  investigation  of 
organic  compound  physisorption  on  model  surface  systems 
of  graphene  and  M0S2.  Effects  of  inclusion  of  a  Si02  sub¬ 
strate  to  model  a  realistic  material  were  also  examined  for 
selected  cases.  Various  adsorbates,  including  those  used  in 
the  experimental  work  of  Perkins  et  i.e.,  electron-donor 
and  acceptor  molecules,  as  well  as  solvents  of  varying  polar¬ 
ity,  were  considered.  Physisorption  characteristics  were  ana¬ 
lyzed  in  detail,  where  optimized  hybrid  material  system 
structures  were  elucidated  for  the  molecular  adsorbates  trie¬ 
thylamine,  acetone  (ACE),  tetrahydrofuran,  methanol, 
2,4,6-trinitrotoluene  (TNT),  o-nitrotoluene,  o-dichloroben- 
zene,  and  1,5-dicholoropentane.  Charge  transfer  analyses  in 
the  optimized  surface — chemical  adsorbate  systems  provided 
partial  rationalization  of  observed  sensing  trends  for  some 
molecules  upon  inclusion  of  a  substrate,  and  have  shown  that 
a  MoS2-based  FET  could  potentially  consist  of  a  more  sensi¬ 
tive  platform  for  electron-donating  chemical  groups.  Indeed, 
local  density  of  states  (LDOS)  results  demonstrated  small 
effects  upon  adsorption,  however,  showing  a  relatively  stron¬ 
ger  response  for  an  electron-donor  upon  inclusion  of  the 
underlying  Si02  substrate. 

II.  COMPUTATIONAL  DETAILS 

DFT  calculations**’*^  were  carried  out  by  solving  the 
Kohn-Sham  equations*^  using  the  Vienna  ab  initio  simula¬ 
tion  package  (VASP).***’^*  The  electron-ion  interactions  for 
Mo,  S,  C,  O,  N,  Cl,  and  H  were  described  by  projector  aug¬ 
mented  plane  wave  (PAW)  potentials. A  plane-wave 
energy  cutoff  of  400  eV  was  used  for  all  calculations,  found 
to  be  sufficient  for  these  systems.**  The  generalized  gradient 
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approximation  (GGA)  exchange-correlation  functional  of 
Perdew-Burke-Ernzerhof  (PBE)  was  used,^"*  and  the  vdW- 
DE^^  non-empirical  correction  was  applied  to  correct  for 
London  dispersion  upon  molecular  adsorption.  Considering 
the  large  size  of  some  of  the  molecules,  a  large  8x8  super¬ 
cell  was  modeled  with  planar  dimensions  of  19.7  x  17.1 
for  SL-graphene  and  25.1  x  21.7  A^  for  SL-M0S2  surfaces. 
Brillouin  zone  sampling  of  the  total  energy  was  based  on  the 
technique  by  Monkhorst  and  Pack^®  with  a  k-point  mesh  of 
3  X  3  X  1 ,  increased  to  6  x  6  x  1  for  LOGS  and  charge  trans¬ 
fer  calculations.  A  k-point  mesh  of  5  x  5  x  1  was  used  for 
structural  optimization  and  total  energy  calculations  for  the 
M0S2  and  graphene  surfaces  on  Si02  with  a  smaller  cell, 
increased  to  10  x  10  x  1  for  LOGS  and  charge  calculations. 
A  k-point  mesh  of  3  x  3  x  1  for  2,4,6-trinitrotoluene  on 
M0S2  and  graphene  surfaces  on  SiG2  was  used  because  a 
larger  cell  was  required.  In  calculation  of  the  total  energy  of 
the  system  for  a  relaxed  configuration,  all  atoms  in  a  super¬ 
cell  were  allowed  to  move  in  all  directions  and  the  structure 
was  relaxed  until  the  forces  acting  on  each  atom  were  con¬ 
verged  to  less  than  0.01  eV/A.  Bader  charge  analyses^^  were 
used  to  analyze  charge  transfer. 

III.  RESULTS  AND  DISCUSSION 

A.  Physisorption  on  free-standing  graphene  and  M0S2 
surfaces 

Eirst,  we  describe  details  of  adsorption  on  graphene  and 
M0S2  surfaces.  The  graphene  surface  structure  is  shown  in 
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Eig.  1(a),  having  C-C  bond  lengths  of  1.42  A.  M0S2  has  a 
closely  packed  layered  structure  with  each  Mo  atom  coordi¬ 
nated  by  six  sulfur  atoms  in  a  trigonal  prismatic  unit,  with 
Mo-S  bond  lengths  of  2.40  A  (Eig.  1(b),  side  views  Eigs.  1(c) 
and  1(d)).  Mo  atoms  in  a  plane  between  two  planes  are  in 
2-H  stacking,^^  with  S  at  a  perpendicular  distance  of  3.16  A. 
Eor  each  of  the  molecular  adsorbates,  namely,  triethylamine, 
acetone,  tetrahydrofuran  methanol,  2,4,6-trinitrotoluene, 
o-nitrotoluene,  o-dichlorobenzene,  and  1,5-dicholoropentane, 
various  adsorption  sites  were  examined,  as  described  below. 
The  adsorption  energy  Eads  was  derived  by  subtracting  the 
total  energy  of  any  given  molecule  E„,oi  and  that  of  the  corre¬ 
sponding  fully  relaxed  unadsorbed  M0S2  or  graphene  surface 
system  E.,,,,/  from  the  total  energy  of  the  hybrid  material  sys¬ 
tem,  namely,  {Emoljsuif^-  ^ads  —  ^mol/suif  ^surf  ^mol- 

Because  the  PBE  functional  does  not  account  for  London 
dispersion  upon  physisorption,  the  reported  geometries  and 
adsorption  energies  were  calculated  using  the  vdW-DE^^ 
non-empirical  correction  to  the  GGA  functional,  previously 
shown,  for  example,  to  be  important  for  peptide  aromatic  res¬ 
idue  adsorption.'®  Results  using  PBE  with  the  vdW-DE  cor¬ 
rection  indeed  demonstrate  a  stronger  relative  adsorption, 
e.g.,  for  2,4,6-trinitrotoluene,  as  expected,  as  summarized  in 
Table  I.  The  optimized  structures  are  depicted  in  Eig.  2. 
Distances  of  heavy  atoms  in  the  molecular  moiety  to  the 
adsorbing  surface  are  listed  in  Table  II.  Because  the  chemical 
structures  are  physisorbed  on  the  surfaces  we  note  that  in  all 
cases  the  internal  structure  of  the  molecule  has  been  hardly 
perturbed. 


FIG.  1.  The  minimum  energy  struc¬ 
tures  of  (a)  graphene,  (b)— (d)  M0S2 
with  two  side  views  (c)  and  (d). 
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TABLE  I.  Adsoiption  energies  (eV)  of  molecules  on  preferred  adsorption 
sites  (structures  shown  in  Fig.  2),  calculated  with  PBE  and  the  vdW-DF  non- 
empirical  correction  to  the  PBE  functional. 


Graphene 

M0S2 

PBE 

PBE, 

vdW-con‘ected 

PBE, 

PBE  vdW-corrected 

Triethylamine 

-0.05 

-0.34 

-0.08 

-0.40 

Acetone 

-0.04 

-0.26 

-0.02 

-0.26 

T  etrahy  drofuran 

-0.05 

-0.30 

-0.03 

-0.31 

Methanol 

-0.04 

-0.14 

-0.04 

-0.16 

2,4,6— trinitrotoluene 

-0.10 

-0.59 

-0.05 

-0.57 

o— nitrotoluene 

-0.06 

-0.42 

-0.04 

-0.40 

o-dichlorobenzene 

-0.04 

-0.37 

-0.02 

-0.37 

1 ,5-dicholoropentane 

-0.01 

-0.46 

-0.02 

-0.48 

Two  sites  for  adsorption  for  triethylamine  on  graphene 
were  analyzed,  with  the  nitrogen  atom  in  both  starting  con¬ 
figurations  on  top  of  graphene’s  carbon  atom,  but  in  one  case 
the  nitrogen  and  ethyl  groups  bent  inwards  and  in  the  other 
outwards.  The  inward  structure  was  slightly  lower  in  energy, 
with  a  binding  energy  of  — 0.34eV  (Table  I)  as  compared  to 
— 0.30eV  for  outward  adsorption.  The  molecule  was  dis¬ 
placed  from  its  starting  position,  with  the  nitrogen  atom 
moving  from  the  (0,  0)  top  site  to  a  position  in  the  xy-plane 
(—0.34,  0.23  A)  (see  structure  in  Fig.  2(a)).  On  M0S2,  the 
site  with  the  nitrogen  atom  and  ethyl  groups  of  triethylamine 
pointing  outwards  was  preferred  by  0.03  eV.  The  molecule 
undergoes  signihcant  motion  and  was  transferred  from  an 
initial  position  on  top  of  a  sulfur  atom  to  a  hollow  site  (Fig. 
2(b)).  Optimized  distances  of  the  nitrogen  atom  to  the  surfa¬ 
ces  are  somewhat  shorter  for  graphene  as  compared  to  M0S2 
(4.44  A  vs.  4.63  A,  respectively,  in  Table  II).  Stronger 
adsorption  is  noted  on  the  M0S2  surface,  of  — 0.40eV. 

In  adsorption  of  acetone,  investigated  for  five  nonequiva¬ 
lent  sites,  a  different  behavior  for  the  two  surfaces  is  shown. 
Adsorption  sites  considered  included  the  oxygen  atom  posi¬ 
tioned  on  the  top,  bridge  and  hollow  sites  and  central  CO  par¬ 
allel  to  the  surface,  and  in  another  case  the  oxygen  atom 
pointed  away  from  the  surface.  The  preferable  adsorption  site 
on  the  graphene  surface  (so-called  CO-flat,  shown  in  Fig. 
2(a))  resulted  in  a  binding  energy  of  — 0.26eV  (Table  I), 
while  a  top  site  with  adsorption  energy  of  — 0.19eV  was  the 
next  strongest  adsorbed  structure.  In  this  case,  the  CO  bond  is 
no  longer  parallel  to  the  surface,  but  the  O  atom  is  slightly 
tilted  toward  the  graphene  C  atom  at  a  distance  of  5.96  A 
(Table  II).  Similarly,  the  CO  bonds  of  acetone  on  the  top, 
bridge  and  hollow  sites  do  not  stay  perpendicular  to  the  gra¬ 
phene  surface  upon  relaxation  but  are  tilted  with  a  slightly 
smaller  angle  than  the  most  preferred  site  (CO-flat).  On 
M0S2,  the  on-top  site  is  slightly  preferred  by  0.03  eV,  as  com¬ 
pared  to  the  CO-flat  site.  Examination  of  the  structures  indi¬ 
cated  small  differences  between  the  two  sites,  with  an  O-S 
bond  distance  (do-s)  of  0.14  A,  which  is  smaller  than  for  the 
most  preferred  adsorption  site.  Fig.  2(b)  shows  the  minimum 
energy  structure,  and  the  corresponding  distance  (Table  II) 
demonstrates  closer  proximity  to  the  surface  than  for 
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FIG.  2.  Top  and  side  views  of  adsorbed  triethylamine,  acetone,  tetrahydrofu- 
ran,  methanol  (left  to  right  top  panel)  and  2,4,6-trinitrotoluene,  o- 
nitrotoluene,  o-dichlorobenzene  and  1,5-dicholoropentane  (bottom  panel  left 
to  right,  respectively)  on  (a)  graphene  and  (b)  M0S2  surfaces. 

graphene.  The  adsorption  energies  for  the  remaining  three 
sites  were  in  the  range  of  —0.13 — 0.16eV. 

For  tetrahydrofuran,  with  out-of-plane  oxygen  atom  in 
the  ring  structure,  we  considered  four  distinct  sites  for 
adsorption.  Two  sites  consider  adsorption  in  a  perpendicular 
direction  through  the  oxygen  atom  on  top  and  hollow  sites 
on  the  surface,  while  the  remaining  two  sites  consider  planar 
adsorption  parallel  to  the  surface  with  the  center  of  the  tetra¬ 
hydrofuran  ring  aligned  with  the  graphene  ring  (planar-top 
site)  and  the  center  of  the  tetrahydrofuran  on  top  of  the  C 
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TABLE  II.  Distances  (A)  of  adsorbed  molecules  for  prefeiTed  adsorption  sites  on  graphene  and  M0S2  surfaces. 


Triethylamine 

Acetone 

Tetrahydrofuran 

Methanol 

2,4,6-trinitrotoluene 

o-nitrotoluene 

o-dichlorobenzene 

1 ,5-dicholoropentane 

Graphene 

dc-ci 

4.89 

3.78 

4.05 

4.21 

4.04 

4.09 

4.02 

3.40 

dG-C2 

4.78 

dc-cs 

4.30 

dc-o 

5.96 

3.41 

3.37 

dc-c 

1.53 

I.5I 

1.54 

1.41 

1.41 

1.40 

1.52 

dc-o 

1.23 

1.43 

1.43 

dc-H 

I. to 

1. 10 

I. II 

I. II 

1.09 

1.09 

1.09 

1.11 

dc-ci 

1.73 

1.80 

dc-ci 

3.40 

dc-N 

1.47 

1.47 

1.47 

^N-O 

1.24 

1.24 

do-N 

4.44 

4.04 

do-H 

0.97 

M0S2 

ds-ci 

4.17 

4.08 

3.92 

4.55 

3.75 

4.06 

3.94 

4.00 

ds-C2 

4.52 

4.01 

ds-C3 

4.08 

ds-o 

3.53 

3.54 

3.31 

dc-c 

1.53 

I.5I 

1.55 

1.42 

1.40 

1.40 

1.52 

dc-o 

1.23 

1.43 

1.43 

dc-H 

I. II 

1. 10 

I. II 

1.11 

1.09 

1.09 

1.09 

1.11 

dc-ci 

1.74 

1.80 

ds-ci 

4.09 

dc-N 

1.47 

1.48 

1.48 

^N-O 

1.24 

1.24 

ds-N 

4.63 

3.78 

do-H 

0.97 

atom  of  graphene  (planar-hollow  site).  Planar  adsorption 
sites  were  of  lower  energy,  and  among  those  the  planar- 
hollow  site  was  slightly  preferred.  The  adsorption  energy  of 
tetrahydrofuran  on  graphene  was  — 0.30eV  (Table  I)  and  the 
structure  is  shown  in  Fig.  2(a).  For  M0S2,  a  similar  structure 
was  found  to  be  a  minimum  energy  structure  (Fig.  2(b))  with 
an  adsorption  energy  of  —0.31  eV  (Table  I).  Distances  to  the 
surfaces  are  generally  comparable  (Table  II),  although  some¬ 
what  shorter  to  the  M0S2  surface. 

Six  adsorption  sites  for  methanol  on  graphene  and  M0S2 
were  investigated,  through  either  the  oxygen  or  carbon  atom 
on  top,  bridge  and  hollow  sites.  In  Fig.  2(a),  we  show  the 
adsorption  of  methanol  at  the  oxygen  atom  on  a  graphene 
hollow  site,  although  adsorption  on  all  studied  sites  was  sim¬ 
ilar,  with  binding  energies  of  ca.  —0.12 — 0.14eV.  However, 
for  the  M0S2  surface,  preference  for  methanol  adsorption  on 
a  hollow  site  through  an  oxygen  atom  was  calculated  (Figure 
2(b)),  with  a  binding  energy  of  — 0.16eV.  Earlier  calcula¬ 
tions  using  a  vdW-DF  correction  to  PBE  showed  comparable 
magnitudes  and  site  preferences  for  adsorption.'^ 

Adsorption  of  2,4,6-trinitrotoluene  and  o-nitrotoluene 
was  modeled  such  that  the  benzene  rings  were  placed  on  top 
of  a  graphene  ring.  During  relaxation,  the  2,4,6-trinitrotolu¬ 
ene  molecule  moved  slightly  inward  as  compared  to 
o-nitrotoluene,  and  a  fully  relaxed  molecule  was  at  a  distance 
of  4.04  A  and  3.78  A  from  the  graphene  and  M0S2  surfaces. 


respectively  (Table  II  and  Fig.  2).  o-nitrotoluene  dislocated 
slightly  in  the  xy-plane  and  moved  outward  in  the  z-direction 
during  relaxation.  The  optimized  structure  was  at  a  distance 
of  4.09  A  from  the  graphene  surface,  with  a  similar  value  for 
M0S2  (Table  II).  Unlike  for  graphene,  the  C-NO2  and  C-CH3 
bonds  did  not  stay  aligned  to  the  Mo-S  bonds  of  M0S2 
(Fig.  2(b)).  Larger  adsorption  energies  were  calculated  for 
stacked  configurations  on  the  surfaces  (Table  I). 

o-dichlorobenzene  adsorption  was  similarly  modeled 
on  graphene  and  M0S2,  initially  placed  at  3.00  A  away  from 
the  graphene  surface.  It  moved  farther  from  the  surface 
(4.02^.05  A)  upon  relaxation,  with  a  binding  energy  of 
—0.37  eV.  For  adsorption  on  single-wall  carbon  nanotubes,  top 
ring  adsorption  was  predicted  to  be  the  most  preferred  site 
with  binding  energies  of  the  same  order. Similar  to  graphene, 
o-dichlorobenzene  prefers  a  top  ring  site  on  the  M0S2  surface. 
The  M0S2  Mo-S  ring  is  larger  than  benzene,  and  therefore, 
o-dichlorobenzene  lies  within  larger  Mo-S  ring  with  a  binding 
energy  similar  to  that  for  graphene  (Table  II  and  Fig.  2). 

Three  sites  for  1,5-dicholoropentane  adsorption  were 
studied  on  graphene,  specifically  adsorbed  parallel  to  the  gra¬ 
phene  surface  with  two,  three  and  five  out  of  the  five  carbon 
atoms  pointing  inward.  We  found  that  the  site  that  has  all 
five  carbon  atoms  of  1,5-dicholoropentane  aligned  parallel  to 
the  graphene  surface  showing  the  strongest  binding  energy, 
of  —0.46  eV,  next  to  the  site  with  three  carbon  atoms  bound 
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TABLE  III.  Accumulation  (+)  and  depletion  (— )  of  charge  (e)  for  molecular 
adsorbates  on  free-standing  graphene  and  M0S2  surfaces. 


Graphene 

M0S2 

Triethylamine 

-0.003  (-0.128)“ 

-0.001  (-0.250)“ 

Acetone 

0.005  (0.037)“ 

-0.039  (-0.147)“ 

T  etrahy  drofuran 

0.002 

-0.021 

Methanol 

-0.001 

0.002 

2,4,6— trinitrotoluene 

0.140(0.017)“ 

0.038(0.031)“ 

o-nitrotoluene 

0.006 

0.009 

o-dichlorobenzene 

-0.001 

0.037 

1 ,5-dicholoropentane 

0.005 

0.012 

“Results  for  Si02-supported  graphene  and  M0S2  surfaces. 


to  the  surface  (— 0.43eV).  This  implies  a  stronger  binding 
energy  with  a  larger  number  of  alignment  of  the  bonds.  We 
found  this  to  hold  for  the  M0S2  surface  too.  A  slightly  stron¬ 
ger  bound  1,5-dicholoropentane  (by  ^0.02  eV)  stayed  rela¬ 
tively  far  from  the  M0S2  surface  with  a  S-C  bond  length 
(ds-c)  of  4.00  A,  as  compared  to  a  graphene-C  distance 
(dc-c)  of3.40A  (Table  II). 

B.  Charge  transfer  upon  chemical  adsorption 

Comparing  changes  in  charge  transfer  upon  molecular 
adsorption  to  the  relative  conductivity  change  of  a  nano-FET 
is  complicated  because  of  consideration  of  only  one  molecule 
in  the  calculations,  and  lack  of  consideration  of  impurities  or 
defects. However,  in  this  work,  we  seek  an  intrinsic 
response  of  the  surface  to  molecular  adsorption.  Bader  charge 
analyses  are  summarized  in  Table  III.  Charge  transfer  on  free¬ 
standing  surfaces  demonstrate  small  changes,  of  0.00-0.04  e 
except  for  2,4,6-trinitrotoluene,  consistent  with  physisorbed 
adsorbates  in  other  work.*”*  Note  that,  for  example,  for  trie- 
thylamine  and  acetone,  these  small  values  are  consistent  with 
the  corresponding  distances  of  the  nitrogen  and  oxygen  atoms 
to  the  surfaces,  respectively  (Table  II).  Triethylamine  and  ace¬ 
tone  had  no  noticeable  effect  on  shifting  the  Dirac  point  away 
from  the  Fermi  level  in  the  LDOS  for  the  graphene  surface 
because  of  the  minimal  charge  transfer  (Fig.  3).  Results  for 
the  M0S2  surface  for  triethylamine  and  acetone  similarly 
show  a  small  response  (Figs.  4(a)  and  4(b)).  On  the  other 
hand,  charge  transfer  to  2,4,6-trinitrotoluene  was  calculated  to 
be  somewhat  larger  than  for  other  adsorbates  (Table  III),  as 
expected.  The  removal  of  charge  from  the  graphene  surface 
upon  adsorption  of  the  strong  electron-acceptor  will  conse¬ 
quently  shift  the  Dirac  point  in  LDOS  away  from  the  Fermi 
level  towards  higher  energies  (conduction  band),  as  shown  in 
Fig.  3,  while  the  LDOS  of  M0S2  in  Fig.  4  also  show  a  small 
shift  around  the  Fermi  level. 

Experimentally,  however,  robust  sensitivity  was  indi¬ 
cated  for  triethylamine,  a  laboratory-safe  decomposition 
product  from  the  V-series  of  nerve  gas  agents,  and  also  ace¬ 
tone,  particularly  for  MoS2.*^  To  explore  potential  interac¬ 
tions  with  the  underlying  Si02  substrate  that  may  explain 
this  behavior,  charge  transfer  for  surfaces  on  silica  was  con¬ 
sidered  for  selected  cases.  A  model  based  on  a-quartz  was 
considered,  because  it  represents  the  local  character  of 
Si02’s  amorphous  structure, as  previously  rationalized.^'  It 
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E-E,(eV) 

FIG.  3.  LDOS  of  carbon  atom  for  free-standing  graphene  (black  line)  com¬ 
pared  with  carbon  atom  of  graphene  with  adsorbed  TEA,  ACE,  and  TNT. 
Inset  shows  the  enlarged  version  of  Dirac  point  shift  for  supported  graphene. 

was  shown  that  differences  in  adsorption  of  graphene  on  a 
Si02  substrate  based  on  a-quartz  or  cristobalite  models  are 
small, and  because  our  goal  was  to  understand  relative 
changes  upon  molecular  adsorption,  only  one  type  of  surface 
was  evaluated.  It  was  pointed  out^*  that  experimental  and 
theoretical  studies  showed  that  the  upper  subsurfaces  of 
cleaved  Si02  undergo  surface  reconstruction  that  results  in 


E  -  E,  (eV) 


FIG.  4.  LDOS  of  (a)  Mo  atom  and  (b)  S  atom  of  free-standing  M0S2  surface 
(black  line)  compared  with  Mo  atom  of  MoS2-suiTace  with  an  adsorbed 
TNT  molecule.  Insets  show  the  enlarged  version  of  the  plot  around  Fermi 
energy  showing  a  shift  in  the  LDOS. 
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oxygen-terminated  surfaces. It  was  also  shown  that  the 
optimized  stmctures  of  graphene  on  oxygen-terminated  or 
hydroxylated  (0001)  Si02  surfaces  are  similar.^'  We  note 


that  comparing  to  an  experimental  /i-type  device/^  presum¬ 
ably  due  to  defects  and  impurities  as  shown  in  recent  theoret¬ 
ical  work,*^  it  is  unlikely  that  there  are  silanol  groups  on  the 
Si02  surface.  We  therefore  considered  an  oxygen-terminated 
a-quartz  Si02  surface. 

The  choice  to  model  the  substrate  by  a  (0001)  a-quartz 
surface  was  additionally  based  on  the  small  lattice  mismatch 
for  both  M0S2  and  graphene  surfaces.  The  optimized  structure 
of  the  a-quartz  Si02  model  structure  resulted  in  Si-0  bond 
lengths  between  1.62-1.65  A  with  a  Si-O-Si  bond  angle  of 
142.6°,  in  agreement  with  earlier  calculations. The  calcu¬ 
lated  lattice  constant  for  graphene  was  2.46  A,  and  therefore, 
a  4  X  4  cell  of  graphene  (32  C  atoms)  was  adsorbed  on  a  2  x  2 
cell  of  (0001)  Si02,  leading  to  a  lattice  mismatch  of  ^  3%.  A 
2x2  cell  with  four  layers  of  a-quartz  Si02  terminated  with 
oxygen  was  modeled  with  a  vacuum  region  of  20  A.  The  lat¬ 
tice  parameter  of  a-quartz  Si02  was  calculated  as  4.85  A  and 
for  M0S2  was  3.20  A.  A  3  X  3  cell  of  M0S2  (9  Mo  and  18  S 
atoms)  was  therefore  adsorbed  on  a  2  x  2  cell  of  (0001)  Si02 
(16  Si  and  40  O  atoms),  leading  to  a  lattice  mismatch 
of~2%.  The  small  strain  was  shown  not  to  affect  the  elec¬ 
tronic  properties.'^  For  computational  feasibility,  the  cell  size 
was  reduced  by  almost  a  half  from  what  we  used  for  the 
unsupported  graphene  and  M0S2  surfaces.  Optimized  distan¬ 
ces  between  the  top  layer  of  Si02  and  the  adsorbed  surfaces 
were  2.85  A  for  graphene  and  2.97  A  for  M0S2  (shown  in  Fig. 
5  for  M0S2).  Charge  densities  for  graphene  and  M0S2  on 
Si02  are  shown  in  Figs.  6(a)  and  6(b),  respectively. 

The  Si02  substrate  causes  p-doping  by  1.071  e  and  0.738 
e  of  graphene  and  M0S2  surfaces,  respectively.  As  a  conse¬ 
quence  of  the  graphene  p-doping,  which  is  consistent  with 
results  for  cristobalite  Si02,*"'  the  Dirac  point  of  the  LDOS  is 
shifted  from  valence  band  towards  the  higher  energies  (con¬ 
duction  band)  by  1.15  eV,  as  shown  in  Fig.  7.  The  LDOS  of 
Mo  and  S  atoms  for  the  M0S2  surface  on  Si02  are  shown  in 
Figs.  8(a)  and  8(b),  demonstrating  a  shift  of  the  band  gap  of 


(b) 

0  0  0 

*.0’  *0 

i\)Hr  w 

0  0  0 

*  * 

•  • 

4r  4 

•  • 

FIG.  6.  Surface  contour  plots  of  charge 
densities  for  graphene  and  M0S2 
adsorbed  on  Si02  plotted  along  a  plane 
passing  through  the  center  of  (a)  car¬ 
bon  of  graphene,  silicon,  and  oxygen 
atoms  and  (b)  molybdenum,  sulfur,  sil¬ 
icon,  and  oxygen  atoms. 
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E-E,{eV) 

FIG.  7.  LDOS  of  carbon  atom  of  free-standing  graphene  (dashed  black  line) 
and  supported  graphene  (solid  black  line)  with  TEA  and  ACE.  Inset  shows 
the  enlarged  version  of  Dirac  point  shift. 

about  0.25  eV,  consistent  with  the  charge  transfer.  This 
response  affects,  in  turn,  changes  in  charge  transfer  upon  mo¬ 
lecular  adsorption,  for  example,  the  electron-accepting 
2,4,6-trinitrotoluene  molecule  shows  accumulation  of  charge 
on  free-standing  surfaces  of  graphene  and  M0S2.  Charge 
transfer,  however,  to  graphene  and  M0S2  surface  with  the 
Si02  substrate  is  smaller,  because  of  the  p-doping  of  the 
surfaces,  as  expected. 


FIG.  8.  LDOS  of  Mo  atom  (a)  and  sulfur  atom  (b)  for  free-standing  M0S2 
(dashed  line)  and  supported  M0S2  with  adsorbed  TEA  and  ACE.  Inset  shows 
the  enlarged  version  of  LDOS  around  Fermi-level. 


On  the  other  hand,  depletion  of  charge  on  triethylamine 
for  adsorption  on  both  graphene  and  M0S2  surfaces  modeled 
on  a  substrate  is  noted,  of  0.128  e  and  0.250  e,  respectively 
(see  Table  III),  as  well  as  for  acetone  on  M0S2  (e.g.,  0.147  e 
on  M0S2,  Table  III,  however,  inconclusive  for  graphene). 
Interestingly,  the  intrinsic  effect  is  stronger  for  M0S2,  quali¬ 
tatively  consistent  with  experimental  observations.^^  A 
smaller  Dirac  point  shift  because  of  a  smaller  charge  transfer 
for  the  system  with  the  triethylamine  and  acetone  on  gra¬ 
phene  is  noted  in  Fig.  7.  Because  triethylamine  and  acetone 
«-dope  the  Si02-supported  M0S2  surface  more  significantly, 
the  LDOS  shift  towards  the  valence  band  (lower  energies) 
upon  triethylamine  adsorption  (Fig.  8)  is  more  noticeable 
compared  to  graphene. 

IV.  CONCLUSIONS 

In  this  work,  we  investigated  by  DFT  calculations  the 
physisorption  of  the  molecular  adsorbates  triethylamine, 
acetone,  tetrahydrofuran,  methanol,  2,4,6-trinitrotoluene, 
o-nitrotoluene,  o-dichlorobenzene,  and  1,5-dicholoropentane 
on  surfaces  of  single-layer  M0S2  and  graphene,  in  order  to 
understand  parameters  that  may  affect  the  intrinsic  response 
for  chemical  sensing  application.  Detailed  investigation  of 
adsorption  sites  for  the  hybrid  material  systems,  and  correction 
for  London  dispersion  in  the  DFT  functional,  assured  accuracy 
of  the  calculated  adsorption  energies.  Charge  transfer  analyses 
and  LDOS  results  for  the  optimized  surface  with  chemical 
adsorbates  demonstrated  very  small  effects.  However,  inclu¬ 
sion  of  an  underlying  Si02  substrate  model  surfaces  in  some 
cases,  provided  rationalization  of  the  observed  sensing  for  trie¬ 
thylamine,  for  example.  These  results  may  assist  in  devising 
chemical  sensors  for  improved  sensitivity. 
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